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Abstract 

The quark production in classical color fields is investigated with a focus on the 
induction of an electromagnetic current by produced quarks. We show that the color 
SU(2) and the SU(3) theories lead significantly different results for the electromag¬ 
netic current. In uniform SU(2) color fields, the net electromagnetic current is not 
generated, while in SU(3) color fields the net current is induced depending on the color 
direction of background fields. Also the numerical study of the quark production in 
inhomogeneous color fields is done. Motivated by gauge field configurations provided 
by the color glass condensate framework, we introduce an ensemble of randomly dis¬ 
tributed color electric fluxtubes. The spectrum of photons emitted from the quarks 
by a classical process is shown. 


1 Introduction 

Photons are one of the most important observables to explore the properties of high-energy 
matter produced in relativistic heavy-ion collisions, since they are penetrating probes 
which are free from strong interactions. In recent years, direct photon measurement by 
the PHENIX [H [2] and the ALICE Collaborations [3l 0] posed a puzzle to theoretical 
understanding of photon production. The direct photon puzzle consists of two issues, 
(i) State-of-the-art theoretical calculations laHlElIH] underestimate direct photon excess 
measured in experiments mm- (ii) A large elliptic flow of direct photons comparable 
in magnitude to that of pions are observed I3I1!. It is difficult to explain these two 
issues simultaneously because a large direct photon yield favors photon emission at early 
time when temperature of the system is high, while a large elliptic flow indicates that 
photons are emitted at later time when the flow of matter is largely developed. This 
puzzle naturally drives people to think about a novel mechanism of photon production, 
e.g. piiniiiiiiiaiis]. 

In Ref. [13], it has been shown that the photon spectra measured in p-|-p, d-|-Au, 
Au-|-Au collisions at RHIC and Pb-|-Pb collisions at LHC satisfy geometric scaling. This 
result indicates that the photons are emitted in early time when the saturation momentum 
Qs is the only dominant scale of the system. In the leading order description of the 
color glass condensate, the system right after a heavy-ion collision consists of longitudinal 
color electric and magnetic fields dSlIHlE]. This state is called Glasma. The photon 
production in the Glasma has been studied in Refs. [giinj. Although promising results are 
presented, these studies are based on a simplified model for the gluon distribution function 
in the Glasma. To obtain deeper understanding of the photon production in the Glasma, 
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calculations based on a first principle is needed. The purpose of the present study is to 
take a first step toward this direction. 

In order to investigate the photon production in the Glasma (the central rapidity 
region of nucleus-nucleus collisions), we need to first compute quark production because 
the initial states described by the color glass condensate (see e.g. [T8l [191 ISHi) purely 
consist of gluons, which do not directly couple to photons. This situation is in contrast 
to proton-nucleus collisions, where the proton can be viewed as high-energy partons, not 
as classical gauge fields. The photon and dilepton production in proton-nucleus collisions 
have been investigated by calculating the processes that a quark in the proton emits a 
photon due to the interaction with strong color classical fields of the nucleus |H1 (221123122] • 

In the leading order description of the color glass condensate and the Glasma, gauge 
fields are treated as classical fields because the coupling is small, g' <C 1, and the fields 
are strong, A ~ l/g( (see e.g. [25] )• The quark dynamics in the leading order is described 
by retarded solutions of the Dirac equation under the classical gauge field. For generic 
inhomogeneous gauge fields, we need to resort to lattice-discretized numerical calculations 
to solve the Dirac equation. In studies of real-time dynamics of fermions on lattice, 
there have been considerable advances [261 Ea ESI EH Eq]. Given these advances, we 
will give a first-principle description of the quark dynamics under classical gauge fields, 
and investigate the induction of electromagnetic (EM) current accompanying the quark 
production. In this paper, we assume a system in a fixed box, namely boost-invariant 
expansion is not taken into account. Besides we focus on only the one-point expectation 
value of the EM current operator. From the expectation of the current, the spectrum of 
photons produced by a classical process can be computed. 

In the context of real-time lattice simulations of the early-time dynamics in heavy-ion 
collisions, the color SU(2) theory is often used instead of the color SU(3) just for simplicity, 
and it is sometimes expected that SU(3)c results would be similar to SU(2)c. Indeed, it has 
shown that the dynamics of classical gauge fields exhibits no qualitative difference between 
SU(2)c and SU(3)c [M]. We will show that, however, this is not the case for the induction 
of EM currents. In SU(2)c, there is a tendency that two color components contribute to 
the EM current with opposite signs to each other. We will show that this cancellation 
is exact in uniform fields. On the other hand, in SU(3)c, the degree of the cancellation 
between different color components depends on the color direction of background fields. 
Thus, the EM current can be several orders larger than that in SU(2)c. 

The rest of this paper is organized as follows. The induction of EM currents by 
uniform color electric fields is discussed in Sec. We compare lattice numerical results 
with known analytic solutions to show the validity of our numerical method. In Sec. 
we investigate the quark production and the subsequent induction of EM currents in 
inhomogeneous background fields, whose configuration mimics the Glasma. The photon 
spectrum associated with a classical emission process is shown at the last. Section is 
devoted to summary and discussions. The lattice numerical method we employ is explained 
in the Appendix. In particular, our method to treat doubler modes is described in detail. 

2 Uniform color electric fields 

In this section, we investigate the induction of EM currents in uniform color electric 
fields. Such a simple configuration for color fields enables us to do analytical studies 
on real-time dynamics of quark production via the Schwinger mechanism [321 I33j . and 
thus clarifies the difference between SU(2) and SU(3) color fields for the induction of 
EM currents. Therefore, studying the uniform color field configurations can provide us 
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with useful information to understand more complicated and realistic situation, which is 
discussed in the next section. 

Furthermore, the uniform field configurations can be used as a benchmark test for the 
real-time lattice numerical computations. By comparing numerical results with known 
analytic results, we can check the validity of our numerical method on lattice before 
applying it to field configurations that analytical solutions are not available. See Appendix 
m for details of the lattice numerical method. 

In this paper, we consider only one family of quark flavor which has an EM charge 
e and mass m, for simplicity. We will briefly comment on the effects of multi-flavors in 
Sec. m 

2.1 Abelianization 

We briefly review a method to simplify the color structure of the Dirac equation under 
uniform color electromagnetic fields |33l I34j . Thanks to this simplification, the induction 
of EM currents in uniform color fields can be understood intuitively. 

If the field strength is spatially uniform, its color dependence can be factorized as 

= ( 1 ) 

where n“ is a constant vector in color space. This kind of fields are called covariant 
constant fields [35l [36]. This field strength is given by the gauge field 

= A^{x)n^ , ( 2 ) 

where A^{x) is related with by F^^ = d^A^i — dyA^. The gauge field Q is a solution of 
the Yang-Mills equation. Because the color vector n“ is constant, the structure of the co¬ 
variant derivative = d^ — igA^{x)n°'T°‘ can be reduced from SU(Ac) to [U(I)]'^‘'. Since 
n“r“ is an W x N^. hermitian matrix, one can diagonalize it by a unitary transformation 
as 

UrFT^U^ = diag (u;i, • • • , . (3) 

The color-rotated quark field V'i = iUtp)^ (i = 1, • • • , W) obeys 

[n'" {9)1 - imgAf,) - m] = 0, (4) 

which has the same structure as the U(l) theory but a coupling e is replaced by the 
effective coupling Wig. 

The SU(2) color space is isotropic since the rank of SU(2) is one. For any unit color 
vector n“, the color matrix n“T“ takes the same form after diagonalization: 

t/n“™t = diagQ,-0. (5) 

Therefore, the effective couplings are always +gl2 and —gj^. 

For Nc> 3, the color space is anisotropic in a sense that the effective couplings depend 
on the direction of the color vector n“ in a gauge-invariant way. For SU(3)c, there are two 
independent diagonal generators. If we take the generators T" = with the Gell-Mann 
matrix A“, the diagonal generators are and T®, by which the diagonalized color matrix 
n“r“ can be expanded: 

UrFT^U"^ = T® cos 0 - T® sin 9. (6) 
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The angle 0 is related with the color vector n“ through the second Casimir invariant 


C2 = [d 




n n n as 


sin"^ 39 = 3 C 2 . 


(7) 


By this equation, the angle 9 characterizes the color direction of the vector n“ in a gauge- 
invariant way. The eigenvalues of the matrix n“r“ are parametrized by 9 as 
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Unlike the SU(2)c case, these effective couplings have the color direction dependence. 


2.2 EM currents in uniform color electric fields 

Having discussed the diagonalization (Abelianization) of the color matrix, one can easily 
prove the cancellation of EM currents in SU(2)c uniform electric fields (see Fig. [^. Let 
us suppose that the electric field is applied to the z-direction with the strength Eq (> 0). 
In the case of SU(2)c, the effective couplings are independent of the color direction of the 
gauge field. A quark with color 1 couples to the electric field with the effective coupling 
of +\g, and it also has an EM charge +e (we assume e > 0). This quark is accelerated 
to the + 2 :-direction by the electric field and causes a positive EM current. An anti-quark 
is accelerated to the opposite direction but it has a negative EM charge, so that it also 
induces a positive EM current. On the other hand, a quark with color 2 has the effective 
coupling —^g. Therefore, this quark is accelerated to —z-direction. Although the color 
charges are opposite in sign between color 1 and 2, the EM charge is independent of the 
color components. Consequently, the quark with color 2 induces a negative current, which 
cancels the current induced by the quark with color 1. The same cancellation happens for 
anti-quarks. In this way, we can conclude that uniform SU(2) color electric fields do not 
induce an EM current. 

Clearly, this argument does not apply to the SU(3)c case. Although there is alway a 
tendency of the cancellation among different color components, the cancellation may not 
be perfect. When 9 = 0 (modulo vr/S), SU(3)c is reduced to SU(2)c because one of three 
quarks decouples from the color field. For other values of the color angle 9, a nonzero EM 
current is induced as we will show explicitly in the next subsection. 
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Figure 1: SU(2)c quarks in a constant color electric field. 
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Figure 2: The time dependence of the EM current induced by the uniform electric field for 
various values of the color angle 9. The lattice numerical results (bold lines) are compared 
with the simple estimation (14) (thin black lines). Parameters used for these numerical 


computations are mjQs = 0.1, Qs^x = Qs^-y = 30, 
N. = 192. 


= 15, Nx = Ny = 48, and 


2.3 Numerical results 


In the following we consider a uniform color electric field turned on at t = 0: 


r?a _ 

-^03 — 


Eon“ 

0 


(t>0) 

{t < 0). 


(9) 


This held conhguration is a solution of the Yang-Mills equation coupled to a pulse-like 
external current 

= -Eon“Y3<5(t). (10) 

In Fig. the expectation value of the EM currents obtained by lattice numerical 
calculations are plotted for some different values of the color angle 9 in the case of SU(3)c. 
Only the z-component of the current is shown because the x- and y-components are zero 
by the symmetry of the system. When 9 = 0, the EM current is indeed zero because 
SU(3)c is reduced to SU(2)c. For other values of 9, a nonzero EM current is generated 
and its magnitude largely depends on the value of 9. This strong dependence of the EM 
current on the color angle 9 stands in contrast to the weak dependence on 9 that the the 
total quark number density shows |37l l33] . This result indicates that EM observables may 
be more sensitive to color held conhgurations than hadronic observables. 

The linear increase of the EM current seen in Fig. and its dependence on 9 can 
be understood from analytical results. Thanks to the Abelianization, the Dirac equation 
under the uniform electric held ([^ can be analytically solved, and the distribution function 
of produced quarks is available in a closed form |33j. Although the distribution function 
has a complicated expression in terms of special functions, its behavior is well described 
by the following simple function [32]: 

/ 77^2 _|_ _j_ p2 \ 

fs,cit;p) = exp (-TT— ^ j 0(p^;0,WcgEot) , (11) 
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where the function 0(x; a, b) is defined by 


0(x; a, b) 


1 min(a, b) < x < max(a, b) 
0 otherwise. 


( 12 ) 


The expectation value of the EM current can be expressed by the analytic expressions for 
the mode functions. However, it still contains the integration over momentum, which needs 
to be evaluated numerically, making the parametric dependence on the color angle hardly 
accessible. Instead of doing the numerical integration, we use the simple expression © to 
evaluate the EM current, and compare it with the results of lattice numerical calculations. 
In terms of the distribution function, the EM current can be expressed as 

4„(t) = 2ej:/^£i/,,(t;p). (13) 

(2-Kyujp 


This is called conduction current, which is caused by movement of charged particles. In 
addition to the conduction current, there is another contribution to the current, called 
polarization current, which cannot be expressed by the normal distribution function. Ex¬ 
cept the times right after the turning-on of the electric field, the polarization current is 
negligible compared with the conduction current under a strong field, gE ^ m? [32] . We 
therefore consider only the conduction current. If we substitute Eq. © into Eq. ( [I^ , 
and assume Pzjojp ~ sgn(p 2 ), which is valid for \wcg\EQt 3> m, we get a simple expression 
for the longitudinal EM current. 




4e 


(2vr) 


'"\'^<^ii\^owcg\wcg\Elt. 
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(14) 


This simple estimation is plotted in Fig. as thin black lines with nonzero constants added 
to fit the numerical data. Eq. (14) well reproduces the 0-dependence that the numerical 


results show, especially at later times. 

In Fig.[^ the coefficient of the linear increase of the current (14) is plotted as a function 
of the color angle 9 for various quark mass. As indicated by Eqs. ([^, this is a periodic 
function with a period of 27r/3. 
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Before closing this section, let us emphasize the fact that the knowledge of the Abelian- 
ization has not been used in the numerical computations. We have just solved Eqs. (A3), 
(A4) and (AlO). Therefore, this numerical calculation can be easily extended to more 
complicated gauge field configurations as we will discuss in the next section. 


3 Glasma-like electric fields 


In the previous section, we have dealt with the uniform electric fields to clarify the physics 
of the quark production in classical non-Abelian gauge fields and the difference between 
SU(2) and SU(3) color fields for the induction of the EM current. In this section, we study 
EM currents induced by more realistic configurations for the background gauge fields. 

The gauge fields in the Glasma are uniform in the longitudinal direction reflecting the 
boost invariance of the system, and are inhomogeneous in the transverse plane |17) . The 
inhomogeneity in the transverse plane is characterized by the saturation scale Qs- One 
of the characteristic features of the Glasma is the existence of the longitudinal magnetic 
fields. However, the existence of magnetic fields at the initial time is inconsistent with 
the assumption of the free initial condition for spinors. In contrast, initial inhomogeneous 
electric fields can be given by zero gauge field A(0) = 0, and thus they can be consistent 
with the free initial condition for spinors. In the T-rj coordinates, one can construct the 
initial condition for the spinors at r = 0 that is consistent with the evolution of the gauge 
fields describing a collision of color glass condensates IMI. In this case, the existence of 
the initial magnetic fields does not cause the inconsistency. However, we cannot employ 
these initial conditions since we assume a system in a non-expanding box. Instead of it, we 
employ the gauge fields that initially consist of only electric fields. To mimic the Glasma, 
we introduce an ensemble of electric fluxtubes. In the transverse x-y plane, we randomly 
distribute Wube = QsLxLy-numhevs of fluxtubes which has a Gaussian profile with mean 
radius 1/Qs and which are uniform in the longitudinal direction: 


Q 


2 Mube 


= 0, ®_l) = — exp n“ , 
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(15) 




where is transverse coordinates of a center of a fluxtube, which is randomly distributed 
in the transverse plane, and n“ is a unit color vector whose orientation is random. The 
distance \x — is measured in a way respecting the periodic boundary condition. The 
electric field ([Tm trivially satisfies the Gauss’s law (A8). 


Figure displays the transverse profiles of the initial energy density of a randomly 
generated electric field (15) for SU(2)c and SU(3)c, respectively. Energy density times 
is plotted since this combination is independent of the value of the coupling g. In the 
following, all dimensionful quantities are shown in a dimensionless combination scaled by 
powers of Qs- In Fig. the time evolution of the gauge field energy density averaged over 
volume is exhibited. Owing to the inhomogeneity in the transverse plane, magnetic fields 
are quickly developed, and the equipartition of energy between electric and magnetic 
components is realized both for SU(2)c and SU(3)c case. As observed in Ref. [3l], the 
SU(2)c and the SU(3)c theories show qualitatively similar behavior for the non-Abelian 
gauge fields. 


Under these classical gauge fields, we have solved the Dirac equation (AlO) and com¬ 


puted the EM current (All). Because the color directions of the background gauge fields 


are random in the transverse plane, the space average of the EM current is nearly zero. 
In Fig. we plot the space average of the square of the EM current for SU(2)c (Left) and 
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Figure 4: Transverse profiles of the initial energy density (times g^) of a randomly 
generated electric field, IQ\. (Density plots with color bars in the right side of each 

figure.) Left; SU(2)c, Right: SU(3)c. Lattice parameters are Qs^x = Qs^y = 10, and 
N^ = Ny = 48. 


SU(3)c (Right) as a function of time. Unlike the uniform field case, a nonzero current is in¬ 
duced even for SU(2)c. However, the current for SU(2)c is several orders smaller than that 
for SU(3)c. In particular, the longitudinal component is smaller than the transverse ones 
in the case of SU(2)c although the initial background electric field is along the longitudinal 
direction. While the argument of the cancellation of the EM currents between two colors 
is strict only for a uniform system, this result indicates that there is still tendency of the 
cancellation even in the inhomogeneous system. The nonzero but relatively small current 
originates in the violation of the Abelianization argument in nonuniform fields. Hence, 
the transverse components are larger than the longitudinal one for SU(2)c. In contrast, 
for SU(3)c, the cancellation is not perfect even in uniform fields, as shown in the previous 
section. Therefore, there is no reason the EM current is suppressed like the SU(2)c case. 
Indeed, for SU(3)c, the longitudinal current is dominant over the transverse ones. Figure 
shows the transverse coordinate dependence of the absolute value of the EM current for 

SU(2), SU(3)e 




Q gt Q gt 

Figure 5: Time dependence of the electric and magnetic components of the gauge field 
energy density averaged over space, Left: 

SU(2)c, Right: SU(3)c. Lattice parameters are the same as Fig. [3 
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Figure 6: Plots of H^QV) a function of time. Left: SU(2)c, Right: 

SU(3)c. Parameters used for these computations are mjQs = 0.1, QsLx = QsLy = 10, 
QsLz = 15, and Nx = Ny = Nz = 48. 
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Figure 7: The transverse coordinate dependence of the absolute value of the EM current 
divided by the EM coupling, | /(cQs), for SU(3)c at the time Qgt = 5. (Density plots 

with color bars in the right side of each figure.) Parameters are the same as those in Fig.|^ 
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Figure 8: The photon energy spectrum for SU(3)c at the time Qgt = 5. The px- 
distribution with hxed py = 0 and the py-distribution with fixed Px = 0 are compared. 
Parameters used for this computation are mfQs = 0.1, QsLx = QsLy = 20, QgL^ = 15, 
Nx = Ny = 64, and = 48. 


SU(3)c at the time Qst = 5. Reflecting the inhomogeneity of the background gauge field, 
the current has inhomogeneity whose typical scale is the order of l/Qs- 

Once an EM current is induced, it can emit photons. Photons emitted from a coherent 
EM current (nonzero one-point expectation of the EM current operator) form coherent 
classical EM helds. This process is simply described by the classical Maxwell equation 
coupled to the EM current; 

— V X B = —Jem • (16) 

We have numerically solved the lattice version of the Maxwell equation with the EM 
current induced by the SU(3)c random fields. Erom the classical EM fields, we can compute 
the photon energy spectrum by the following equation: 


with 


dEpfi 

(PpTdz 


1 

2(24)2 


E{Pt) + B(p-p) 


(17) 


E(Pt) — 

j(f XT , 

(18) 

B(Pt) = 

[cfxTB{xT)e-^^T-=^^ - 

(19) 


Because the system is uniform in the z-direction, we have not done the Eourier transform 
in z. If we do it, the p^-distribution becomes the delta function S{pz). Eigure shows 
the photon energy spectrum 0 at the time Qgt = 5. As the background color field is 
characterized by the scale Qs , the photon energy spectrum has a width determined by the 
same scale Qs- In the figure, pa;-distribution with fixed py = 0 and py-distribution with 
hxed Px = 0 are compared. Rehecting the inhomogeneous structure of the background, 
the two distributions differ especially at the low momentum region pt ^ Qs- 
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Pt/Qs 


Figure 9: The photon number spectrum (20) as a function of pT- Results of two runs 
with different random initial backgrounds and the average over four runs are plotted. 
Parameters used for these computations are the same as those in Fig. A function 
3.2 X 10“^ g-2.73pT/Qs obtained by a fitting at the region 1 < < 4 is plotted as well. 


In Fig. the photon number spectrum defined by 

d^ph _ 1 dFipfe 

(Ppxdz oj dPpxdz ’ 

with OJ = pt = -^Px + Py^ i® plotted as a function of px- The pT-distribution is computed 
by taking average over points {px,Py) which give the same px- To show the dependence 
on the random background field configurations, the results of two runs with different 
initial color field configurations are plotted, as well as the average over four runs. At the 
low momentum region, px < Qs, the spectrum fluctuates relatively largely run by run, 
in accordance with the observation we made in Fig. At the high momentum region, 
the dependence on the random filed conhgurations is relatively small, and the spectrum 
shows an exponential behavior. In Fig. a function proportional to ^ which is 

obtained by a fitting at the region 1 < < 4, is also plotted. It is worth noting that 

this exponential function resembles a thermal distribution although the source of 

the photons is far away from thermal equilibrium. 

4 Summary and discussions 

In this paper, we have investigated the induction of electromagnetic (EM) currents as¬ 
sociated with the nonperturbative quark production from classical color fields. In weak 
coupling and strong field regimes, which may be realized in the initial stage of heavy-ion 
collisions, the quark production in the leading order with respect to the coupling g is 
described by the retarded solutions (mode functions) of the Dirac equation under the clas¬ 
sical gauge field. Once the mode functions are obtained, the EM current can be computed 
from them. We have numerically solved the Dirac equation on real-time lattice. 

To clarify the significant difference between the SU(2)c and the SU(3)c theories for 
the induction of EM currents, we have considered, at first, uniform and constant color 
electric fields as background fields. This field conhguration enables us to simplify the 
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Figure 10: The leading order diagrams for the photon 2-point function {A^{x)Atj{x)). 
Solid liens denote quark propagator dressed by background non-Abelian gauge fields. Left: 
Disconnect diagram that represents the classical process of photon emission, which we have 
studied in this paper. One loop tadpole corresponds to the one point expectation of the EM 
current operator, (JemC^^))- Right: Connected diagram that describes a genuine quantum 
process of photon production, which is not considered in this paper. 


color structure of the theory by diagonalization. Thanks to this simplification, it can be 
easily shown that a net EM current is not generated in uniform SU(2) color fields since the 
contributions from two colors cancel each other. In contrast, in uniform SU(3) color fields 
a nonzero net EM current is induced depending on the color direction of the background 
electric field. 

After confirming that our lattice numerical calculation in the uniform color field repro¬ 
duces known analytic solutions, we have applied our numerical method to more realistic 
inhomogeneous backgrounds. Motivated by the Glasma, we have introduced fluxtube-like 
color electric fields which are randomly distributed in the transverse plane [Eq. (15)]. Un¬ 
der such inhomogeneous gauge fields, the argument of the exact cancellation of the EM 
current for SU(2)c cannot be applied, and thus nonzero EM current is generated even 
for the SU(2)c case. However, the EM current for SU(2)c is several orders smaller than 
that of SU(3)c, because there is still a tendency of the cancellation between two colors. 
Therefore, it is important to pay attention to the difference between SU(2)c and SU(3)c, 
when one studies the induction of EM current and subsequent photon production. Using 
the EM current obtained by numerical computations, we further calculated the EM fields 
induced through the Maxwell equation. Erom the EM fields, the photon spectrum has 
been computed. Because the background color fields are characterized only by a single 
scale Qs-, the photon spectrum has a width given by the same scale Qs- Reflecting the 
random nature of the background color fields, the photon spectrum fluctuates relatively 
largely at small momentum region, px < Qs- We have found that the photon spectrum 
shows, at the region pt^ Qs, an exponential behavior with — Qs/3. 

There are several directions of possible extensions of this study. 

While we have dealt with fixed geometry, the systems of relativistic heavy-ion collisions 
show nearly boost-invariant expansion in the longitudinal beam direction. The boost- 
invariantly expanding system can be conveniently described by the r-rj coordinates. Once 
we take account of the expanding geometry, we can utilize the McLerran-Venugopalan 
model [381 139] for the initial condition of the classical background fields instead of the 
random fields we introduced by hand. The initial condition for the quark mode functions 
which is compatible with the background fields provided by the color glass condensate 
framework has been derived in Ref. m- 

It is also interesting and important to go beyond the leading order description in 
the coupling g. The next-to-leading order contributions for the quark production can 
be computed by adding fluctuations to the initial condition for the classical gauge fields 
and considering the back reaction from produced quarks to the background color fields 
[29] . If the gauge field fluctuations are involved, the background gauge field is no longer 
uniform in any direction. In this case, the numerical cost to solve the Dirac equation 
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becomes unfavorably expensive, and thus we may need to resort to the statistical method 
to sample the mode functions. 

In this paper, we have considered only one family of quark flavor. Let us comment on 
how our results are modified by the inclusion of multi-flavors in the case of the uniform 
electric field. In the leading order of the strong field and weak coupling approximation, 
quarks with different flavors are produced independently. Therefore, we just need to sum 
up the EM currents contributed from each flavor. The quark masses shown in Fig. 
m/yJgEQ = 0.05,0.65,2.1, correspond to strange, charm, bottom, respectively, if we as¬ 
sume the field strength is related with the saturation scale as qEq = , and take the value 

Qs = 2 GeV. Because the line of m/^/gEo = 0.05 is overlapped by that of m/^/gEo = 0, 
strange and lighter flavors have the same yield. Since the sum of the EM charges of up 
(-|-2e/3), down (—e/3), and strange (—e/3) is zero, the total EM current induced by these 
three light flavors is vanishing. As a consequence, there remains only the contribution 
from charm quark. Since this argument is based on Eq. (14) which is valid only for the 
uniform electric field, consideration of a realistic situation is needed to make a robust 
conclusion. 

The most important extension of this study would be the computation of a genuine 
quantum contribution for photon production. The photon spectrum can be related with 
the photon two point function {A^{x)Ai,{y)). In the leading order with respect to the 
coupling e, there are two kinds of diagrams, shown in Fig. 10, In this paper, we have 


computed the classical EM fields by solving the Maxwell equation which couples to the 
one-point expectation of the EM current operator, and the photon spectra have been 
calculated from the classical EM fields. This classical contribution corresponds to the 
disconnected diagram shown in the left of Fig. 10 The connected diagram shown in 


the right represents a genuine quantum process in a sense that it is related with the 
connected two-point current correlation function {J^yi{x)J'^yi{y))c-, and it is impossible to 
construct it from the classical current {JEMix))- Although the connected contribution 
is the same order in the coupling e as the disconnected one, it is not computed in this 
paper. It is important to investigate this genuine quantum contribution for the photon 
production, in addition to the classical contribution studied in this paper. In particular, 
the cancellation between different colors/flavors may not happen in this contribution, 
because the two vertexes are connected in the right diagram. However, the computation 
of the genuine quantum contribution in the real-time non-equilibrium quantum field theory 
involves several problems, which demand further investigation. The present paper provides 
a basis for these future investigations. 
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Appendix A Numerical method 

A.l Lattice discretization 

Except a few examples, analytic solutions of the Dirac equation under inhomogeneous 
gauge fields are hardly available. Therefore, one needs to resort to numerical methods. In 
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order to numerically solve the classical Yang-Mills equation and the Dirac equation, we 
discretize the 3-dimensional space V = x Ly x into x Ny x A^j;-lattice sites. The 
space coordinates are labeled by 


X — {cLxTlxi CLyTly, ClzTlz) , 


(Al) 


where n* = 1, 2, • • • , Yj and a* = Li/Ni (i = x, y, z). We assume the periodic boundary 
condition. The plane wave is discretized as 




= exp 



/ kxTix I kzTLz \ I 


(A2) 


where ki = -|- 1, • • • , —1,0,1, • • • 

We take the temporal gauge Ag = 0 and keep the time as a continuum variable. In this 
case, the gauge fields are described by link variables Ui{x) = exp {igaiAi{x)) and electric 
fields Ei{x). The equations of motion for these variables are 


dtUi{x) = igaiEi{x)Ui{x) , 


(A3) 


and 


1 1 1 

dtEi{x) = -^ {Im [Ui,j{x) + Ui-j{x)] - (trace)} 

OYd 

where Uij{x) and Ui-j{x) are the plaquettes defined by 

Uij{x) = Ui{x)Uj{x + i)U}{x + j)U'^-{x), 
Ui,.j{x) = U^ix)U]{x + i- j)Uj{x - j)Uj{x - j ), 


and “—(trace)” means 


1 


A - (trace) = A - ^trA • 1 (for A E SU(Ac)). 

’ r. 


The lattice Yang-Mills equations (A3) and (A4) preserve the Gauss’s law 

3 


— Ei{x) - u} (x - i)Ei{x - i)Ui{x - i) 


i=l 


di - 


= 0 . 


(A4) 


(A5) 

(A6) 

(A7) 


(A8) 


Since derivative appears linearly in the Dirac equation, the derivative for the Dirac 
field must be replaced by the center difference to maintain the hermiticity of the theory: 


1 


Diil){x) = ^ Ui{x)'4^{x + i) - U}(x - i)tp(x - i) 

Zdi L 


(A9) 


Using this covariant difference, we need to solve the Dirac equation for the mode functions, 

[ij^do + ifDi - m] = 0 , (AlO) 


as well as the lattice Yang-Mills equations (A3) and (A4). We assume the null initial 
condition for the gauge field: Ui{to) = 1. Hence the quark mode functions must obey the 
free initial condition. 
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On the lattice, the vacuum expectation value of the EM current operator J*(x) = 
e'iijj{x)Yiijj{x) can be expressed in terms of the mode functions as 


«^EM e ^ Y ^ R-® 

s,c p 


'ipp,s,cixh'’Ui{x)i^p^,^^{x + i) 


(All) 


Note that only the negative energy mode functions are necessary to compute such expec¬ 
tations. 

The numerical cost to compute the mode functions V'p,s,c(^) is proportional to {N^NyNzf' 
because the mode functions depend on both x and p. To reduce this rather expensive 
numerical cost, one can employ a Monte Carlo method |27ll30j . However, in this paper, we 
use the deterministic method because we treat gauge fields which are uniform in one direc¬ 
tion at least, of which nature reduces the numerical cost. Under a gauge field uniform in 
the 2 ;-direction, the ^-dependence of the mode functions V’p,s,c(®) i® always There¬ 

fore, the ^-dependence need not be computed numerically. The scaling with the lattice 
size for the numerical cost to compute the mode functions is modified to {NxNy)‘^Nz. A 
merit of the deterministic method is that it is free from statistical errors. This is beneficial 
especially when one computes local expectations under inhomogeneous backgrounds. 


A.2 Treatment of doublers 


As is well known, naively discretized fermionic fields suffer from the problem of doublers. 
The lattice momentum associated with the center difference 


1 

Pi = — sin 
CLi 



(A12) 


has 2-fold degeneracy. It has two zeros at fej = 0 and ki = Ni/2 as Fig. 11 shows. In 
3-dimensional momentum space, there is 8-fold degeneracy in total. We call the region 


■^phys — ^ (fex) ky, kz) 


Ni 

ki — -;-h 1 


Ni 

0 ... — 


(A13) 
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the physical region, and the remain the doubler region Ajoubier- A way to decouple the 
doubler modes from the dynamics is adding the Wilson term to the Dirac equation: 


(i7°5o + * 7 *A - m) ijj{x) + ^ X/ ~ Ui{x)'il^{x + i) - 2il;{x) + Uj (x - i)ip{x - i) 


i=l 


= 0 . 
(A14) 


However, we do not employ the Wilson term by the following reason. 

Because the Wilson term contains the gauge link, it contributes to the EM current. In 


addition to the normal contribution (All), we have 


J* = — 

'^EM,W ^ Y 




(A15) 


In the continuum limit, the contribution from the Wilson term vanishes since ^{x)'il){x) 
is real. However, with a nonzero lattice spacing and under a nontrivial gauge field con¬ 
figuration, this contribution stays nonzero and contaminates the normal contribution. In 
particular, for the EM current, the normal contribution may be small or zero due to the 
cancellation between different colors as discussed in Sec. and thus the Wilson term con¬ 
tribution can be the same order to the normal contribution. To circumvent this problem, 
we employ an alternative way to decouple the doubler modes. 

Because we directly handle the mode functions with the specific initial conditions, 
we can restrict the momentum p of the mode function to the physical modes, 

p G Aphys- By this, the mode functions which have the momentum index p G Ajoubler 
can be simply excluded from the dynamics. However, the effects of the doubler modes 
may not be completely eliminated by this procedure since the mode functions have also 
x-dependence. If we Eourier-transform the negative-energy mode functions as 

'^P,s,cit^p') = , (A16) 

X 


then the momentum p' spans over the full momentum space Aphyg © Adoubier- There is 
still a possibility that the doubler modes contaminate the mode functions through the 
x-dependence. At the initial time, the mode functions are proportional to the plane wave 
e~iP'^, so that ^|J is diagonal in the momentum space: 


'^p,s,cii = 0) pO ot ^P,P' ■ (A17) 

Because we have restricted the momentum p to the physical modes, also the momentum p' 
is in the physical region. In this sense, the effects of doublers are completely eliminated at 
the initial time just by the limitation p G Aphys- By the time evolution, the contamination 
by the doubler modes may come in. This point becomes apparent if we Eourier-transform 
the Dirac equation: 


(f7°5o -m) %^s,c{t,p') 


+ ‘Et" i E Uj{q) - U ](-,)] p'-q) = 0. 

j = l J q 


(A18) 


where Ui{q) is the Eourier transform of the link variable. If the gauge field is inhomoge¬ 
neous, Ui{q) can carry nonzero momentum and the gauge field may kick out quarks from 
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Figure 12: Plots of Ylk' ^{'t'ikx,ky = 0;k'^,ky) for (a) k^ = 0 and (b) kx = 12. This 
quantity equals to at the initial time Qst = 0. Shaded regions represent the doubler 

region. Parameters used for this computation are rajQs = 0.1, QsLx = QsLy = 10, 
QsLz = 15, and Nx = Ny = Nz = 48. 


the physical momentum region to the doubler region. By explicitly computing 
we can monitor how the contamination occurs. 

For that purpose, we consider the following quantity: 

M{t-,kx,ky-,kx,ky) = ^ (^ 19 ) 

^ fc.eAphy, fc' s,c 

where kx and ky {k'x and k'y) are integers corresponding to the momenta px and py {p'x and 
p'y), respectively. Because 'il^p^s,c{^iP') is always diagonal with respect to pz and p'z under 
longitudinally uniform gauge fields we consider in this paper, the sum over the longitudinal 
momenta is already taken. Thanks to the normalization condition for the mode functions, 
the function M satisfies 

Y Mit-,kx,ky;kx,ky) = 1. (A20) 

Jc' h' 

At the initial time, the function M is totally diagonal, 

M[0] kx, ky', kx;, ky) = 6kx,k'x.^ky,ky ■ (-^ 21 ) 

It therefore vanishes if {k'x.,ky) G Ajoubier since {kx,ky) G Aphys- Figure IT2| exhibits how 
the function M evolves in time under the inhomogeneous gauge field introduced in Sec. 
Both of the strength and the extent of the inhomogeneity of this gauge field is characterized 
by the scale Qs- As examples, Ylk' ^{kikx,ky',k'x.,ky) with {kx,ky) = (0,0) and (12,0) 
are shown. At the initial time Qgt = 0, the fc(,-distribution is localized at k'^. = kx- Since 
the inhomogeneous background gauge field carries nonzero transverse momentum, the /c(,- 
distribution gets broadened by time evolution. However, the extent of the broadening is 
an order of Qs, which is a typical momentum the background field carries, and far smaller 
than the cutoff momentum. In particular, when k'x is near to the border with the doubler 
region (the case oi k'x = 12, plotted in the right), the broadening is weaker than the case k'^. 
is at a low momentum region. As long as we restrict {kx, ky) to the physical region, {k'x, k'y) 
at the doubler region is not significantly excited. In particular, the zero momentum point 
of the doubler {k'x = 24 in this case) is not at all activated. Although we have shown only 
two examples of the parameters {kx, ky), the broadening is weak for any value of {kx, ky), 
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at least within the timescale of the order of l/Qs; which we are interested in. To conclude, 
by restricting the momentum p in the index of the mode functions '4’p,s,c{^) to the physical 
region, we can suppress the contamination by doubler mode^ This would be not true 
if the background field is highly inhomogeneous, or if we continue the time evolution till 
a late time. For example, if one reaches thermal equilibrium after a long time evolution, 
the physical modes and the doubler modes would be equally excited unless one uses other 
method to suppress doublers. However, in the background field and time scale we are 
interested in, the contamination by the doublers does not matter. 

This method has an advantage in numerical cost. Because we restrict the momentum 
p to the physical region, which has the volume of 1 /2^ of the whole momentum space, the 
numerical cost is reduced by factor eight. 
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